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Abstract 

Background: Abscisic acid (ABA) can regulate tine expressions of many stress-responsive genes in plants. However, in 
defense responses to pathogens, mounting evidence suggests that ABA plays variable roles. Little information exists 
about genome-wide gene expression in ABA responses in tomato {Solanum lycopersicum L.), a model fruit crop plant. 

Results: Global transcriptome profiles of tomato leaf responses to exogenous ABA were generated using lllumina 
RNA-sequencing. More than 173 million base pair reads were mapped onto the tomato reference genome and the 
expression pattern differences between treated and control leaves were assessed. In total, 50,616 transcripts were 
generated. Among them, 42,583 were functionally annotated in the NCBI non-redundant database and 47,877 in 
the tomato genome reference. Additionally, 31,107 transcripts were categorized into 57 functional groups based on 
Gene Ontology terms, and 14,371 were assigned to 310 Kyoto Encyclopedia of Genes and Genomes pathways. In 
both the ABA treatment and control samples, 39,671 transcripts were available to analyze their expressions, of 
which 21,712 (54.73%) responded to exogenous ABA. Of these transcripts, 2,787 were significantly differently 
expressed genes (DEGs). Many known and novel ABA-induced and -repressed genes were found. Exogenous ABA 
can influence the ABA signaling pathway with PYR/PYL/RCARs-PP2Cs-SnRK2s as the center. Eighteen PYL genes 
were detected. A large number of genes related to various transcription factors, heat shock proteins, pathogen 
resistance, and the salicylic acid, jasmonic acid, and ethylene signaling pathways were up-regulated by exogenous ABA. 

Conclusions: The results indicated that ABA has the potential to improve pathogen-resistance and abiotic stress 
tolerance in tomato. This study presents the global expression analysis of ABA-regulated transcripts in tomato and 
provides a robust database for investigating the functions of genes induced by ABA. 

Keywords: Tomato, Exogenous ABA, RNA-Seq, ABA signaling pathway, Transcription factors, Heat shock proteins, 
Pathogen-related proteins, ROS scavenging enzymes 



Background 

Tomato {Solanum lycopersicum L.) is an important fruit 
crop and a model system in plants. The basic chromo- 
some number of tomato is 2n = 24, and wild forms range 
from diploids to hexaploids [1,2]. The genome of the to- 
mato was recently published; it possesses about 35,000 
genes, a rich resource to be studied by scientists [3]. 
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The hormone abscisic acid (ABA) regulates numerous 
developmental and functional processes, including stoma- 
tal aperture and hydraulic conductivity, seed dormancy, a 
key phase transitions throughout the plant lifecycle [4-6]. 
An important role of ABA in a variety of plants is to in- 
crease tolerance to stresses such as drought, salinity, cold, 
and heat [7-9]. ABA was recognized as an important sig- 
naling molecule that can trigger signaling and regulating 
mechanisms to cope with adverse stresses. Under drought 
conditions, plants synthesize and redistribute ABA, lead- 
ing to guard cell responses that close stomata, thereby re- 
ducing plant water loss [10]. Research in chickpea {Cicer 
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arietinum) has shown that ABA can improve cold toler- 
ance by, in part, improving the water status of leaves and 
anti-oxidative ability and can enhance heat tolerance via 
the accumulation of osmoprotectants [11,12]. Rice {Oryza 
sativa) seeds pretreated with ABA had improved salt- 
tolerance, with lower Na^ and CI" levels and Na^/K^ ratio 
as well as higher and Ca^* concentrations, proline ac- 
cumulation, soluble sugar content, and grain yield [13]. 

Mounting evidence suggests that ABA plays contrasting 
roles during defense responses to pathogens [14,15]. It in- 
terfered at multiple levels with pathogen stress signaling 
and suppressed or promoted phenylalanine ammonia lyase 
(PAL) activity [16,17]. The different plant species assayed 
in these studies can account, in part, for the various ob- 
served effects of ABA on pathogen resistance. ABA is as- 
sociated predominantly with pathogen susceptibility in 
tomato, but both negative and positive effects were re- 
ported in Arabidopsis [18]. Moreover, the synergistic or 
antagonistic roll of ABA in disease resistance is universally 
acknowledged to depend on the type of pathogen, its 
specific way of entering the host, the timing of the defense 
response, and the type of tissue affected, as well as the 
concentration of ABA [19]. ABA responses during biotic 
and abiotic stresses primarily affect levels of gene ex- 
pression [7] . ABA-induced genes mainly encoded proteins 
associated with dehydrins and enzymes that detoxify react- 
ive oxygen species, enzymes of compatible solute metabol- 
ism, a variety of transporters, regulatory proteins such as 
transcription factors (TPs), protein kinases and phospha- 
tases, and enzymes involved in phospholipid signaling. 
ABA-repressed genes usually encode proteins that are 
associated with growth, including cell wall, ribosomal, 
plasma membrane, and chloroplast proteins [8]. 

To comprehensively understand the functions of ABA 
and related genes, high-throughput screening techniques 
have proven quite effective [20-23]. Matsui et al. (2008) 
completed an Arabidopsis transcriptome analysis after 2 
and 10 h of drought, cold, high-salinity, and ABA-treatment 
conditions using a tiling array, and the results indicated that 
approximately 16% of the genes in the transcriptomic li- 
brary were significantly regulated by ABA. In addition, there 
was greater crosstalk between ABA and drought and high- 
salinity stress signaling processes than between ABA- and 
cold-stress signaling. Similar results were discovered in ex- 
pression profiles of rice genes under various stresses and 
ABA application using a cDNA microarray and RNA gel- 
blot analyses [20]. The impact of ABA on gene expression 
far exceeds that of other plant hormones. About 3,000 
genes detected by microarray were significantly regulated 
by ABA in Arabidopsis seedlings; this number was twice 
and 23 times as many genes as were influenced by methyl 
jasmonate and gibbereUin, respectively [24]. 

With the rapid development of next-generation se- 
quencing technology, RNA deep-sequencing (RNA-seq) 



becomes more efficient and less expensive [25]. The re- 
sults of RNA-seq are highly reproducible, both technic- 
ally and biologically [26,27]. In recent years, numerous 
transcriptome data sets have been produced, made pub- 
licly available, and reanalyzed by other researchers. How- 
ever, so far, few RNA-seq analyses of the genome-wide 
responses of genes to ABA in plants have been reported. 

We investigated ABA responses in tomato, and here re- 
port the results of a comparative transcriptome analysis of 
exogenous ABA-treated tomato leaves using RNA-seq 
technology. The goals were to (i) construct a tomato leaf 
transcriptome; (ii) compare and analyze the transcripts in 
control and ABA-treated plants; and (iii) gain insight into 
stress tolerance and pathogen-resistance induced by ABA 
in tomato. This study presents the transcriptome of to- 
mato leaves responding to ABA and provides a genetic re- 
source that can be used for crop improvement. 

Results and discussion 

Global transcriptome analysis 

Tomato seedling sprayed with 7.58 |imol L''^ ABA solution 
were chosen to study ABA responses. Twenty-four hours 
later, the young third leaves of randomly-selected plants 
from both the ABA-treated and control groups were col- 
lected, labeled as aid (1 d after ABA-treatment) and eld 
(1 d, control), respectively. The tomato ABA/control RNA 
samples were used for deep sequencing on an lUumina 
HiSeq 2000 platform. Each read in the Solexa paired-end 
(PE) sequencing was 101 bp in length. Sequencing gener- 
ated 266.98 million reads, a total raw dataset of 26.96 Gb. 
After trimming, 212.78 million clean reads remained, corre- 
sponding to 20.95 Gb clean data. The dataset of each sam- 
ple, including eld and aid, was represented by over 100 
million reads, a read density sufficient for the subsequent 
quantitative analysis of genes. The sequence reads were 
aligned to the tomato reference genome using SOAPa- 
ligner/soap2 software (http://soap.genomics.org.cn), aUow- 
ing two base mismatches [28]. Of the total reads, 81.97% 
matched either to a unique (36.53%) or to multiple (45.44%) 
genomic locations (Table 1). 

After aligning and assembling, 27,597 genes and 35,051 
transcripts were identified in the transcriptome. Of those 
spliced transcripts, 9,138 (26.07%) matched completely 
with the annotated tomato genome, 12,501 (35.67%) were 
potentially novel isoforms, 3,339 (9.53%) were unknown 
or intergenic transcripts, and 1,634 (4.66%) mapped to the 
complementary strand of an annotated gene. These results 
suggested that some transcripts were probably generated 
from alternative mRNA splicing or were new transcripts. 

The genes and transcripts resulting from Illumina se- 
quencing were merged with the annotated reference gen- 
ome to generate 37,633 genes and 51,606 transcripts using 
Cuffmerge [29]. We removed transcripts with lengths less 
than 150 bp, and the remaining 50,770 transcripts (98.38% 
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Table 1 Number of reads sequenced and mapped to the tomato genome 





cid (control) 


aid (ABA treatment) 


Sum 


Raw bases (bp) 


12,814,237,438 


14,141,031,018 


26,955,268,456 


Raw reads 


1 26,873,638 


140,109,218 


266,982,856 


Clean bases (bp) 


9,936,610,838 


11,016,491,741 


20,953,102,579 


Clear reads 


100455991 


111323699 


211779690 


Total alignment (percent of clean reads) 


82394883 82.02% 


91194594 81.92% 


173589477 81.97% 


Unique matches (percent of clean reads) 


36742770 36.58% 


40622954 36.49% 


77365724 36.53% 


Multi-position match (percent of clean reads) 


45652113 45.44% 


50571640 45.43% 


96223753 45.44% 


Unmatched (percent of clean reads) 


18061108 17.98% 


20129105 18.08% 


38190213 18.03% 



*The number of unique mapping reads plus multimapping reads equals the total number of total alignments. 



of all transcripts), corresponding to 37,093 genes, were 
used for functional annotation and expression profOing in 
each sample. Transcripts of length 150-200 bp accounted 
for 8.14% of the total, those 200-600 bp for 32.53%, those 
600-1000 bp for 20.28%, those 1000-1600 bp for 20.72%, 
those 1600-2200 bp for 9.44%, those 2200-3000 bp for 
5.41%, and those >3000 bp for 12.27%, as shown in 
Figure 1. 

Annotation, functional classification, and KEGG analysis of 
all detected transcripts 

Of the 50,770 transcripts, 45,704 (90.02%) possessed open 
reading frames (ORF), while ORFs in the other 5,066 
(9.98%) were not predicted (Additional file 1: Table SI). Of 
all transcripts, 47,877 were described in the tomato gen- 
ome, and 42,583 had homologs in the NCBI non- 
redundant (NR) protein database. Many of these genes 
were reported to respond to drought, cold, high salinity, 
or ABA, and included both regulatory and functional pro- 
teins. Transcription factors (1,799 transcripts) were the 
largest group of regulatory proteins (Additional file 1: 
Table SI). This result was consistent with the previous 



reports in Ambidopsis and rice using microarray and tOing 
arrays analysis [20,23]. 

For the global functional analysis, all identified tomato 
transcripts were classified into different functional cat- 
egories using Blast2GO (version 2.3.5) [30]. A total of 
31,107 transcripts could be annotated in Gene Ontology 
(GO) and were classified into 57 functional groups, in- 
cluding 23 groups in biological process, 19 in cellular 
component, and 15 in molecular function (Figure 2). 
Within biological process, "metabolic process" (GO: 
0008152) with 17,827 transcripts and "cellular process" 
(GO: 0009987) with 18,169 transcripts were predominant. 
In the category of cellular component, the three main 
groups were "cell" (GO: 0005623, 18,493 transcripts), "or- 
ganelle" (GO: 0043226, 14,423 transcripts) and "cell part" 
(GO: 0044464, 18,493). The categories "binding" (GO: 
0005488) and "catalytic activity" (GO: 0003824) were most 
common in molecular function, represented by 15,938 
and 15,869 transcripts, respectively. 

All detected transcripts were blasted to STRING 9.0 
for further annotation based on Cluster of Orthologous 
Groups (COG) of protein categories. A total of 18,885 
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Figure 1 Length distribution of tlie transcript sequences detected in the ABA-treated tomato leaf transcriptome. 
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biological_process molecular_function 



Figure 2 Gene Ontology (GO) functional annotation of transcripts. All 31,107 transcripts were assigned to at least one GO term and were 
grouped into three main GO categories and 57 groups, 23 groups in biological process, 1 9 in cellular component, and 1 5 in molecular function. 
The right-hand Y-axis represents the number of genes in a sub-category. The left-hand Y-axis indicates the percentage of a specific sub-category 
of genes in each main category. 
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COG categories as listed on the right of the histogram, and the Y-axis indicates the number of transcripts. 
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transcripts with COG annotations were grouped into 25 
functional categories. The largest category was "General 
function prediction only" (3,454 COG annotations, 18.29% 
of 18885), following by "Transcription" (2,021, 10.70%), 
"replication, recombination and repair" (1,776, 9.40%), 
and "signal transduction mechanisms" (1737, 9.20%). In 
addition, only 474 (2.51%) COG annotations belonged 
to the "Function unknown" category (Figure 3). 

KEGG annotation results were retrieved from KEGG data- 
base based on sequence similarity, and 14,371 transcripts 
were assigned to 310 KEGG pathways. The pathways most 
strongly represented were "metabolic pathways" (koOllOO, 
3,476 transcripts), "biosynthesis of secondary metabolites" 
(koOlllO, 1,715), foUowed by "ribosome" (ko03010, 541), 
"plant hormone signal transduction" (ko04075, 457), "starch 
and sucrose metabolism" (ko00500, 267), and "plant-patho- 
gen interaction" (ko04626, 280) (Additional file 1: Table S2). 

Analysis of expressed genes 

In total, 39,671 (78.14% of 50,770) transcripts expressed in 
the tomato leaf transcriptome (Additional fOe 1: Table S3), 
including 38,626 in eld and 37,989 in aid, respectively. 
There were 1,682 transcripts expressed only in the ABA- 
treated group (aid), 1,045 transcripts expressed only in 
the control group (eld), and 36,944 transcripts expressed 
in both libraries, which indicated that ABA may activate 
or repress a fraction of unique transcripts. 

Based on expression level, the transcripts were divided 
into five groups (Table 2). The group with expression 
levels of 1-10 FPKM (Fragments Per Kilobase of exon 
model per Million mapped fragments) was largest, repre- 
senting 44.64% of all transcripts, while those with expres- 
sion levels of <1 FPKM accounted for 28.87%, those with 
10-100 FPKM 23.49%, those with 100-1000 FPKM 
2.73%, and those with >1000 FPKM only 0.26%. Among 
all expressed transcripts, 17,959 showed no changes in ex- 
pression levels (|log2 fold-change (log2FC)| < 0.25). The 
expressions of 21,712 (54.73%) transcripts were altered by 
exogenous ABA (14,559 were up-regulated, 7,153 were 
down-regulated) (Additional file 1: Table S3). Of the tran- 
scripts with altered expression, 2,787 (12.84% of 21,712) 



Table 2 Number of expressed transcripts in control (cid) 
and ABA-treated (aid) libraries 



FPKM 


cid (control) 




aid (ABA treatment) 




Number 


Percent 


Number 


Percent 


>1000 


107 


0.28% 


90 


0.23% 


100-1000 


1117 


2.94% 


975 


2.52% 


10-100 


8964 


23.60% 


9036 


23.39% 


1-10 


16779 


44.17% 


17428 


45.12% 


<1 


11022 


29.01% 


11097 


28.73% 


Total 


37989 




38626 





were significantly changed (|log2FC| > 1 and false discovery 
rate (FDR) < 0.05) and labeled differentially expressed genes 
(DEGs), including 1,952 that were up-regulated and 835 
that were down-regulated (Figure 4). Among the DEGs, 
there were more up-regulated transcripts than down- 
regulated ones, indicating that many genes responded posi- 
tively to ABA treatment. This result was similar to those 
reported in previous studies in Arabidopsis and rice 
[31-33]. 

Molecular functional classification based on GO was also 
performed on the 2,787 DEGs. The up- and down- 
regulated DEGs annotated in GO were grouped into 54 
and 38 groups on G0.1evel3, respectively (Additional file 1: 
Table S4) and into ten classes on G0.1evel2. The most 
common categories were "catalytic activity" (654 up- 
regulated, 295 down-regulated) and "binding" (661 and 
260, respectively), followed by "transmembrane trans- 
porter activity", "substrate-specific transporter activity", 
"transcription factor activity", "translation factor activity", 
and "signal transducer activity". 

Genes related to ABA signaling transduction 

The ABA signaling transduction pathway mainly includes 
four core regulatory components: ABA receptor/pyrabac- 
tin resistance protein 1/PYR-like protein (RCARs/PYRl/ 
PYLs), protein phosphatases type 2C (PP2C), the sucrose 
nonfermentingl-related protein kinase 2 (SnRK2) and 
ABA responsive element binding factors (ABF) [34]. In 
the tomato transcriptome, 18 PYLs, 23 PP2Cs, 12 SnRK2s, 
and 18 ABFs were identified (Table 3, Additional file 1: 
Table S5). 

The discovery of PYR/PYL/RCAR receptors in these 
samples greatly advances the understanding of the ABA 
signaling pathway. Among 18 PYLs in the tomato tran- 
scriptome, the expression level of eight did not differ 
(|log2FC| <0.25) in ABA treatment compared with con- 
trol. The other 10 PYLs showed slight differences in 
their transcript abundance, from 0.55 to 1.43 fold, with 
five up-regulated and five down-regulated genes. Two 
up-regulated PYLs showed moderate transcript level dif- 
ferences (TCONS_00012484, change from 14.16 FPKM 
in cid to 18.07 in aid, TCONS_00036624, from 14.79 to 
17.26 FPKM). The other eight differentially-expressed 
PYLs were found at low abundances (FPKM <10). The 
results indicated that exogenous ABA had a slight im- 
pact on PYL expression because the ABA-mediated sig- 
naling cascade is initiated by the perception of ABA 
receptors [35]. The slight change in upstream ABA- 
receptor gene expression could initiate responses in a 
series of downstream genes. RCARs/PYRl/PYLs were 
recently found to comprise 14 members in Arabidopsis 
thaliana [6,36]. AU PYL family members are ABA-binding 
proteins and that regulate the target PP2Cs in a com- 
binatorial manner [37]. However, we detected 18 PYL 
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Figure 4 Differential expression analysis of all transcripts in the control (cid) and ABA-treatment (aid) libraries. Transcripts tinat differed 
by less than 20% (|log2FC| < 025) were assumed to not change in expression level. Transcripts that satisfied the conditions of "FDR < 0.05" and 
"]log2 fold-change (log2FC)| >1" were considered significantly differentially-expressed genes (DEGs). Other transcripts are noted as "up-regulated 
slightly" or "down-regulated slightly". 



members in this tomato transcriptome, enriching the PYL 
receptor gene family. 

Among the 23 PP2Cs in this study, 13 transcripts chan- 
ged their expression levels (|log2FC| > 0.25), treated by 
ABA, including five that were up-regulated and eight that 



were down-regulated. Of them, only one transcript with 
low levels (TCONS_00014920) was significantly higher than 
the control by 2.05-fold (from 0.37 to 0.74 FPKM). Three 
PP2Cs showed significantiy decreased expression in ABA- 
treatment (aid) by 2.42-fold (TCONS_00014686, from 



Table 3 Genes related to ABA signaling pathway and transcription factors in tomato transcriptome 





Description 


No. up-regulated 


No. down-regulated 


No. unchanged 


Sum 


ABA signaling transduction 


PYP/PYL 


5 


5 


8 


18 




PP2C 


5 


8 


10 


23 




SnRK2 


1 


4 


7 


12 




ABF 


7 


3 


8 


18 


bZIP 


ABRT 


2 


0 


2 


4 




bZIP 


40 


21 


46 


107 


BHLH 


DREB 


2 


5 


4 


11 




BHLH 


47 


22 


66 


135 


MYB 


MYB 


66 


43 


84 


193 




MYC 


6 


2 


11 


19 


AP2 


AP2 


11 


4 


11 


26 




AP2/ERF 


20 


14 


21 


55 




RAV 


1 


0 


2 


3 


NAC 




23 


15 


25 


63 


WRKY 




33 


13 


35 


81 


HSF 




16 


8 


17 


41 


MADS-box 




21 


23 


17 


61 



RCARs/PYRI/PYLs, regulatory components of ABA receptor/pyrabactin resistance proteinl/PYR-Like proteins; PP2Cs, type 2C protein phosphatases; SnRK2, sucrose 
nonfermentingl-related protein l<inase 2; ABF, ABA responsive element binding factors; ABRT, ABA responsive transcription factor; bZIP, basic leucine zipper 
transcription factor; DREB, dehydration responsive element binding protein; bHLH, basic helix-loop-helix transcription factor; MYB, myeloblastosis oncogene; MYC, 
myelocytomatosis oncogene transcription factor; AP2/ERF, apetala 2/ethylene-responsive element-binding factor; NAC, no apical meristem/ATAF/CUP-shaped coty- 
ledons; HSF, heat shock factor. 
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13.79 to 5.70 FPKM), 3.92-foId (TCONS_00012476 from 
16.86 to 4.30 FPKM), and 7.22-fold (TCONS_00030068 
from 2.12 to 0.29 FPKM). The result agreed with the previ- 
ous reports that PP2Cs were negative regulators of ABA 
signaling [8]. In the absence of ABA, PYR/PYL/RCARs are 
not bound to PP2Cs, and PP2C activity is high, which pre- 
vents SnRK2 activation. In the presence of ABA, the com- 
bination of ABA and PYR/PYL/RCARs bind to and inhibit 
PP2Cs, which allows the accumulation of phosphorylated 
SnRK2s and subsequent phosphorylation of ABFs [8]. 

Among the 12 SnRK2 transcripts, five differed slightly in 
level compared with control. Of them, one showed increased 
expression (TCONS_000 11840 from 4.41 to 5.34 FPKM) 
and four showed decreased expression (TCONS_00017720, 
TCONS_00017721, TCONS_00044604, TCONS_00046427). 
Although TCONS_00017721 had moderate transcript abun- 
dance (FPKM from 40.56 to 26.80 FPKM), the other 
three down-regulated SnRK2s were expressed at low levels 
(FPKM < 1.04). SnRK2s are categorized into subclasses I, II 
and III. Subclass I members are not activated by ABA, and 
subclass II and III members are weakly and strongly acti- 
vated by ABA, respectively [38]. Among 10 SnRK2 members 
in Arabidopsis and rice, SAPK8, SAPK9, and SAPKIO were 
in subclass III, and SAPKl to SAPK7 were in subclass I and 
II [38]. In this transcriptome, the only induced SnRK2 
(TCONS_00011840) that was annotated as serine/threonine 
protein kinase SAP/<8-like protein was of subclass III. The 
other four repressed SnRK2s were possibly SAPKl, SAPK3, 
or SAPK7 based on their NCBI blastx matches. Kobayashi 
et al. (2004) demonstrated that only three SnRK2 family 
members were activated by ABA but all were activated by 
hyperosmotic stress, and there were no members activated 
only by ABA. This study also suggested that a small percent- 
age of SnRK2s participated in the signaling transduction 
under exogenous ABA stress. 

In the 18 ABFs, seven were up-regulated and three were 
down-regulated. TCONS_00038922 (from 0.16 to 0.43 
FPKM) and TCONS_00038921 (from 0.43 to 0.10 FPKM) 
were sigruficantiy altered. This result indicated that ABFs 
were mainly positive regulators of ABA response. The mem- 
bers of the ABF gene family belongs to a subfamily of bZIP 
TFs [6]. Most ABFs are known induced in vegetative tissue 
under drought stress and ABA treatment in Arabidopsis [9]. 
However, in the transcriptome, ABFs were expressed at low 
levels in both libraries, so transcription regulator factors 
may target a cascade action [39]. 

Overall, these results agreed with what is known of the 
ABA regulation pathway [6,8,40]. Thus exogenous ABA 
can activate the signaling pathway, and it affected PYLs, 
PP2Cs and SnRK2s, and ABFs. Moreover, members of these 
gene families had different responses to ABA, indicating 
that they probably played positive or negative regulatory 
roles in the ABA pathway. Furthermore, many genes in 
these families had no differences in expression abundance, 



implying that there exists functional redundancy in the 
ABA pathway or that these genes participate in an ABA- 
independent pathway [8]. 

Transcription factors induced by ABA 

With respect to TFs, the results of this study supported 
those of previous research. The TF families of bZIP (111 
transcripts), bHLH (146), MYB (212), AP2/ERF (84), MAC 
(63), and WRKY (81) were identified in this study. Moreover, 
41 Heat shock factors (HSFs) and 61 MADS-boxes were also 
detected (Table 3). In the bZIP family, ABA-responsive TF 
(four transcripts) and bZIP (107 transcripts) were two sub- 
groups. Dehydration-responsive element-binding protein 
{DREB, 11) and bHLH (135) belonged to the bHLH family. 
The MYB family included MYB (193) and MYC (19). The 
AP2 family was divided into three subgroups, including AP2 
(26 transcripts), AP2/ERF domain-containing TF (55), and 
RAVI (3). Reportedly, the most common classes of regula- 
tory sequences induced by ABA are the G-box ABREs 
recognized by members of the bZIP family, and many ABA- 
regulated genes also contain binding sites for proteins of the 
MYB and MYC families. NAC, WRKY, bHLH, and Zn-finger 
classes were also found to participate in some ABA re- 
sponses [8]. 

Although TFs were up-regulated after ABA treatment, 
they had low abundances in both libraries (FPKM < 10) 
(Table 3, Additional file 1: Table S5). The TFs interacted 
with cis-elements in the promoter regions of several re- 
sponsive genes and thus control the expressions of many 
downstream genes, triggering cascade reactions of many 
physiological processes and controlling biochemical reac- 
tions in plant cells [39]. Thus, a slight alteration in the 
transcript abundance of TFs can result in a substantial 
change in downstream gene expression and physiological 
responses. Therefore, although in this tomato transcrip- 
tome ABA affected TFs with low transcript abundances, 
the potential results cannot be neglected. Exceptionally, 
dehydration-responsive element binding proteins {DREBs) 
that bind to drought-response elements (DRE) cis-elements 
were substantially down-regulated (five down-regulated and 
two up-regulated). Reportedly, most DREBs are independ- 
ent of ABA [41]. 

Genes related to lieat sliock protein (IHSP) 

Heat shock proteins (HSPs), including HSP70s, HSP90s, 
HSPlOOs, HSP60s (cpn60s), and small heat-shock proteins 
(sHSPs), are stress-responsive proteins that function as 
molecular chaperones, protecting plants from damage 
under stress [42]. In the tomato transcriptome, 83 tran- 
scripts were identified as HSPs, including 12 HSP90s, 36 
HSPVOs and 24 sHSPs, mainly HSP20s, as well as 11 other 
HSP transcripts (Table 4, Additional file 1: Table S5). 

HSP90s and HSP70s are important members of the 
HSP protein family. The major role of HSP90 is to 
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Table 4 Genes related to heat shock proteins (HSPs), reactive oxygen species (ROS) scavenging, and pathogens 
resistance 

Trait 

Heat shock proteins 



ROS scavenging system 



pathogens resistance 





Description 


No. of up-regulated 


No. of down-regulated 


No. of no change 


Sum 




Hsp90 


3 


1 


8 


12 




HSP70 


8 


8 


20 


36 




sHSP 


4 


11 


9 


24 




Other HSP 


2 


0 


9 


12 




SOD 


1 


1 


8 


10 




CAT 


5 


0 


2 


7 


GSH-AsA cycle 


GLR 


7 


16 


12 


35 




APX 


4 


3 


13 


20 




MDAR 


1 


0 


4 


5 




DHAR 


2 


1 


5 


8 




GR 


1 


0 


3 


4 


GPX pathway 


GST 


23 


20 


33 


76 




GPX 


0 


1 


7 


8 




POD 


33 


28 


31 


92 


PrxR/Trx pathway 


Trx 


25 


29 


79 


133 




PrxR 


1 


5 


4 


10 


Some proteins 


PAL 


6 


3 


7 


16 




PPO 


4 


0 


1 


5 




GLU 


14 


12 


36 


62 




chitinase 


15 


4 


17 


36 


SA signaling 


TGA 


10 


2 


14 


26 




PR1 


15 


12 


18 


45 


JA signaling 


JAR1 


7 


2 


6 


15 




JAZ 


6 


6 


6 


18 


ET signaling 


ETR 


6 


5 


3 


14 




ERF/EREBP 


5 


4 


6 


15 



APX, ascorbate peroxidase; CAT, catalase; DHAR, dehydroascorbate reductase; ERF/EREBP, ethylene response factor; ET, ethylene; ETR, ethylene receptor; GLR, 
glutaredoxin; GLU, beta-1,3-glucanase; GPX, glutathione peroxidase; GR, glutathione reductase; GSH-AsA, glutathione-ascorbate; GST, glutathione-S-transferase; 
HSP, heat shock protein; JA, jasmonic acid; JAR1 , jasmonate resistant 1 ; JAZ, jasmonate ZIM-domain protein; MDAR, monodehydroascorbate reductase; PAL, 
phenylalanine ammonia-lyase; POD, peroxidase; PPO, polyphenol oxidase; PR1, pathogenesis-related proteins; PrxR, peroxiredoxin; ROS, reactive oxygen species; 
SA, salicylic acid; sHSP, small heat shock protein; SOD, superoxide, TGA, TGACG/as-1, binding; Trx, Thioredoxin. 



manage protein folding [43,44], and it also plays key roles 
in signal-transduction networks, cell-cycle control, protein 
degradation, and protein trafficking [45]. HSP70s function 
in preventing aggregation, assisting refolding, protein 
import and translocation, signal transduction, and tran- 
scriptional activation [42,46]. In this transcriptome, three 
up-regulated and one down-regulated HSP90s were 
weakly expressed. Two of the three up-regulated HSP90s 
had high transcript abundances (TCONS_00021180 from 
95.04 to 110.50 FPKM, TCONS_00032910 from 868.34 to 
966.19 FPKM), while the down-regulated HSP90 was 
expressed at low abundance. Among 36 HSP70s, 16 tran- 
scripts differed in expression level, with eight up-regulated 
and eight down-regulated, although the differences in ex- 
pression were slight. The up-regulated TCONS_00037214 
(from 989.52 to 1116.60 FPKM) and the down-regulated 



TCONS_00017664 (from 105.12 to 80.45 FPKM) showed 
high and moderate transcript abundances, respectively. The 
other changed HSP70s were expressed at low abundances. 

Although many HSP70s and HSP90s were slightly up- 
regulated, some showed high transcript abundances. For 
a highly-abundant transcript, a slightly changed ratio 
may substantially alter its abundance and induce physio- 
logical changes. This research indicated that HSP70 and 
HSP90 may play an important role in ABA responses. 

Small HSPs are a large family present ubiquitously in all 
organisms. Plants have many more sHSPs than other eu- 
karyotes [47]. In 24 sHSPs, 15 slightly changed their expres- 
sion levels by ABA, including four up-regulated and 11 
down-regulated transcripts. All up-regulated sHSPs had low 
transcript abundances in both the ABA treatment and con- 
trol. Among the down-regulated sHSPs, there were five with 
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moderate abundance (10 < FPKM < 100) in both libraries 
(TCONS_00036382, TCONS_00026922, TCONS_00028974, 
TCONS_00035725, TCONS_00023776), and the other sfac 
had low abundances. sHSPs bind denatured proteins and 
prevent their irreversible aggregation, so they are referred to 
as 'paramedics of the cell' [48]. Therefore, they have a cru- 
cial and fundamental role in plant cell biology [47,48], and 
contribute to the molecular adaptation to a wide range of 
environmental stresses, including heat, cold, drought, salin- 
ity, and oxidative stress [42]. 

At present, studies of ABA responses to HSP are limited. 
HSPlOl mRNAs in wheat leaves were reported to be in- 
duced by 2 h of dehydration and 50 i^M ABA [49]. In an 
examination of nine rice HSP genes under abiotic stresses 
and ABA treatment, only OsHSP71.1 was induced by ABA 
while OsHSP24.1 was suppressed [50]. In this study, sev- 
eral HSP70s and HSP90s were up-regulated, while most 
sHSPs were down-regulated. These observations imply 
that the HSP genes may play different roles in plant devel- 
opment and abiotic stress responses. Furthermore, the re- 
lationship between ABA and HSPs will be interesting to 
investigate in the future. 

Genes related to the reactive oxygen species (ROS) 
scavenging system 

Overproduction of reactive oxygen species (ROS) and 
overexpression of various antioxidant enzymes in the ROS 
scavenging system occur during almost all biotic and abi- 
otic stresses [51]. In total, there were 408 transcripts iden- 
tified as encoding enzymes in the ROS scavenging system 
(Table 4, Additional file 1: Table S5). They were catego- 
rized into the glutathione-ascorbate (GSH-AsA) cycle (72 
transcripts), the glutathione peroxidase (GPX) pathway 
(176), the peroxiredoxin/thioredoxin (PrxR/Trx) pathway 
(143), superoxide dismutase (SODs, 10) and the catalases 
(CATs, 7). The largest families were genes encoding thior- 
edoxins (Trxs, 133 transcripts), peroxidases (PODs, 92), 
and glutathione S -transferases (GSTs, 76). 

In a stress environment, the first step in O2 reduction 
produces diffusible hydroperoxyl (HO2) and superoxide 
(O2) radicals [52]. Meanwhile, SODs are induced to rapidly 
disproportionate O2 into oxygen and H2O2. However, of 
nine SODs, one with high abundance (TCONS_00025638, 
from 204.01 to 265.26 FPKM) was slightiy down-regulated, 
and another (TCONS_00025639, from 0 to 0.012 FPKM) 
showed slightly increased expression level. The other six 
SODs showed no changes (|log2FC| < 0.25). It was also 
found that the expression levels of SODs showed slightiy 
change after ABA treatment in Arabidopsis [23] . 

CAT reduces H2O2 directly to water and oxygen. In 
this study, the most of the genes encoding CAT were 
up-regulated by ABA, suggesting the importance of 
CATs in scavenging ROS under ABA stress. Among the 
seven identified CAT genes (five up-regulated, two 



unchanged), the abundances of TCONS_00003149, 
TCONS_00003243, TCONS_00049411, TCONS_00051177, 
and TCONS_00003244 in ABA treatment were higher than 
in control by 3.33, 2.20, 1.52, 1.45, and 1.23 fold, respect- 
ively. In particular, TCONS_00049411 (from 115.88 to 
176.01 FPKM) and TCONS_00051177 (from 91.12 to 
131.86 FPKM) had high abundance, and TCONS_00003243 
was expressed a moderate levels (from 9.43 to 20.70 
FPKM). 

Moreover, H2O2 can be also indirectiy scavenged by the 
ascorbate-glutathione (AsA-GSH) cycle, the GPX pathway, 
and the PrxR/Trx pathway [53,54]. The AsA-GSH cycle 
contained 35 glutaredoxins {GLRs, seven up-regulated, 16 
down-regulated), 20 ascorbate peroxidases {APXs, four 
up-regulated, three down-regulated), five monodehydro- 
ascorbate reductases {MDAHs, one up-regulated, four 
unchanged), eight dehydroascorbate reductases {DHARs, 
two up-regulated, one down-regulated), and four gluta- 
thione reductases {GRs, one up-regulated, three un- 
changed). Most transcripts showed slight differences in 
their expression levels between the two libraries. Only 
three transcripts, TCONS_00040877 {GLR, from 3.77 to 
8.24 FPKM), TCONS_00036965 {APX, from 0.00020 to 
1.83 FPKM) and TCONS_00020530 {DHAR, from 0.28 
to 1.44 FPKM) were significantly up-regulated. 

In the GPX pathway, there were 76 GSTs (23 up- 
regulated, 20 down-regulated), eight glutathione peroxi- 
dases (GPXs, one down-regulated, seven unchanged), and 
92 PODs (33 up-regulated, 28 down-regulated). All up- 
regulated genes had low transcript abundances except 
TCONS_00018678 (POD), which was moderately abun- 
dant. Most down-regulated transcripts were expressed 
at low abundance, six occurred at moderate abun- 
dance (3 GSTs: TCONS_00037287, TCONS_00037288, 
and TCONS_00037289; one GPX: TCONS_00034706 
and 3 PODs: TCONS_00018679, TCONS_00009154, and 
TCONS_00004157), and two GSTs (TCONS_00025295 
and TCONS_00039081) were highly abundant. More- 
over, most transcripts in this group showed slight 
alterations in expression level between the libraries. Ex- 
ceptionally, one GST (TCONS_00005473) and three PODs 
(TCONS_00046741, TCONS_00004986, and TCONS_ 
00003494) were significantly up-regulated and three GSTs 
(TCONS_00008673, TCONS_00005786, and TCONS_ 
00025296) and two PODs (TCONS_00003495 and 
TCONS_00009061) were significantly down-regulated. 

There were 133 transcripts encoding Trxs (25 up- 
regulated, 29 down-regulated) and 10 transcripts encod- 
ing PrxRs (one up-regulated, five down-regulated) in the 
PrxR/Trx pathway. All up-regulated Trxs were expressed 
at low abundance in both libraries (FPKM < 10), except 
TCONS_00023279 and TCONS_00024016. Among the 
down-regulated PrxR and Trx genes, more than half 
were moderately or highly abundant. 
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The results of this study showed that most transcripts 
of APX, MDAR, DHAR, GR, GPX, GLR, PrxR, and Trx 
were down-regulated or did not change in expression 
level. These enzymes reduce oxidative damage by H2O2 
in cells directly or indirectly, depending on the GSH 
redox state, which is catalyzed by GR [54-56]. Previous 
research demonstrated that H2O2 and the GSH/GSSG 
couple may interact with other signaling pathways that 
regulate the expression of antioxidant genes during 
stress responses [57]. However, this study showed that 
exogenous ABA repressed the expression of the genes 
encoding enzymes related to the GSH/GSSG couple, 
suggesting that ABA signaling may be independent of 
GSH redox signaling. 

Furthermore, more than half of the GSTs, which may 
catalyze the conjugation reaction of GSH and protect 
SH-containing enzymes against oxidation, were induced 
by ABA [58]. ABA probably also caused the down- 
regulation of genes encoding enzymes that depend on 
GSH. POD catalyzes an early step in lignin synthesis in 
the phenylalanine metabolism pathway, which can in- 
crease the rigidity of cell walls to protect plants from 
pathogenic bacteria [59] and allowing them to resist 
wilting and mechanical weakness imparted by abiotic 
stresses [60]. This study found 33 ABA-induced POD 
genes that maybe participated in these processes. 

Genes related to pathogen resistance 

The phytohormone ABA plays a multifaceted and pivotal 
role in plant immunity. In the whole transcriptome, a large 
number of genes related to disease resistance were de- 
tected, including 16 PAL transcripts (six up-regulated, 
three down-regulated), five polyphenol oxidase (PPO) tran- 
scripts (four up-regulated, one unchanged), 62 beta-1, 
3-glucanase (GLLT) transcripts (14 up-regulated, 12 down- 
regulated), and 36 chitinase transcripts (15 up-regulated 
and four down-regulated). There were more up-regulated 
isoforms than down-regulated ones, and the expression 
changes of up-regulated isoforms were greater than down- 
regulated one. The results suggested that ABA has a 
greater positive than negative impact on tomato pathogen 
resistance. 

The four up-regulated PALs (TCONS_00038804, 
TCONS_00037017, TCONS_00038803, and TCONS_ 
00013292), the up-regulated PPO (TCONS_00036090), 
the up-regulated chitinase (TCONS_00024439), and one 
down-regulated GLU (TCONS_00032191) showed high 
or moderate abundances. Other transcripts of PAL, 
PPO, GLU, and chitinase genes were expressed at low 
abundance, including the significantly up-regulated 
PPO (TCONS_00036092), GLUs (TCONS_00045778 and 
TCONS_00048495), and chitinase (TCONS_00011298) 
(Table 4, Additional file 1: Table S5). PAL, PPO, GLU 
and chitinase are often associated with improved plant 



resistance to pathogens [61-63]. Our previous study showed 
that foliar spraying of exogenous ABA (7.58 \imo\ IS^) 
could enhance resistance to Alternaria solani (early blight) 
infection in tomato by activating defense genes and enhan- 
cing the activities of defense-related enzymes, including 
PAL, PPO and POD, PRl, GLU, and SOD [17]. 

In this transcriptome, many genes in the salicylic acid 
(SA), jasmonic acid (JA), and ethylene (ET) signal trans- 
duction pathways were detected. These pathways play 
well-established roles in phytopathogen defenses [64]. 
TGACG/as-1 bindings (TGAs, 26 transcripts in this 
study) and pathogenesis-related proteins (PRls, 45 tran- 
scripts) are critical component in the SA signaling path- 
way, while jasmonate resistant 1 (JARls, 15 transcripts) 
and jasmonate ZIM-domain proteins (JAZs, 18 tran- 
scripts) are crucial in the JA signaling pathway and 
ethylene receptors (ETRs, 14 transcripts) and ethylene 
response factor/ethylene-responsive TPs (ERF/EREBPs, 
15 transcripts) are vital to the ET signaling pathway 
(Table 4, Additional file 1: Table S5). Most of these genes 
were up-regulated by exogenous ABA and most were 
expressed at low abundance. 

Salicylic acid, a key phytohormone in plant-pathogen in- 
teractions, is typically associated with systemic expression 
of pathogenesis-related protein genes against biotrophs 
and hemibiotrophs. This response often culminates in a 
localized cell death termed the hypersensitive response 
[65-67]. JA and ET predominantly protected Arabidopsis 
against biotrophic pathogens [68]. In tomato, the compara- 
tive transcriptome results indicated that ABA played a 
positive role in biotic stress resulting from the crosstalk be- 
tween ABA and SA, JA, and ET. 

In general, ABA is thought to play a negative role in 
plant resistance to both biotrophic and necrotrophic 
fungi and bacteria, and ABA interacts antagonistically 
with SA, JA, and ET in various ways [69]. For example, 
in tomato, the ABA-deficient mutant sitiens was much 
more resistant to Erwinia chrysanthemi, and exogenous 
ABA (100 |iM) made both the wild type and sitiens more 
susceptible to E. chrysanthemi [18,70]. ABA (100 i^M) 
decreased the disease resistance of rice to rice blast 
{Magnaporthe grisea) by suppressing some genes in the 
SA signaling pathway [71]. However, in this study, ABA 
largely up-regulated the genes related to biotic stresses, 
indicating that ABA has potential to promote pathogen 
resistance in tomato. Furthermore, the concentration of 
ABA applied here (7.58 |iM) was much lower than that 
used in most previous reports (usually 100 \iM). We 
hypothesize that the ABA concentration is very import- 
ant in determining whether plant resistance to patho- 
gens will be improved or suppressed. 

ABA not only mediates tolerance to adverse environ- 
mental conditions, but also affects the resistance to biotic 
stresses. It plays critical roles in the interrelationship 
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between biotic and abiotic stress responses. Further in- 
sights into the core network of ABA signaling in plants 
will be crucial for genetically engineering and breeding 
crop species with both improved abiotic stress tolerance 
and pathogen resistance [72]. ABA has also been recog- 
nized to crosstalk with the "developmental" hormones 
auxin, gibberellic acid, cytokinins, and brassinosteroids in 
plant immunity [14,15,69]. In this study, the expressions 
of many genes related to these hormones were also af- 
fected by ABA (Additional file 1: Table SI). These genes 
may directly or indirectly contribute to plant pathogen 
resistance. 

Quantitative real-time-PCR validation of differentially- 
expressed transcripts from RNA-seq 

To confirm the accuracy and reproducibility of this Illu- 
mina RNA-seq result, 25 transcripts were chosen ran- 
domly for quantitative real-time (qRT) PGR. Those genes 
were involved in metabolism, information transfer, or were 
hypothetical proteins, and included up-regulated, down- 
regulated, and unaffected transcripts. The primer se- 
quences, gene functions, and FPKM and qRT-PCR values 
are listed in Additional file 1: Table S6. The qRT-PCR re- 
sults generally agreed (84%) with the changes in transcript 
abundance determined by RNA-seq, suggesting the reli- 
ability of the RNA-seq data. 

Conclusion 

Genome-wide transcriptome analysis results indicated that 
exogenous ABA can influence the ABA signaling pathway 
with the core of PYR/PYL/RCARs-PP2Cs-SnRK2s. Exogen- 
ous ABA up-regulated many genes related to stress toler- 
ance and pathogen resistance, including various TPs, 
HSP90S, HSP70S, CATs, GSTs, PODs, PALs, PPOs and 
chitinases, as well as the genes TGA, PRl, JARl, JAZ, ERT 
and ERF/EREBP in the SA, JA, and ET signaling pathways. 
These results suggested that ABA has the potential to im- 
prove the abiotic and biotic tolerance of tomato. The study 
extends the knowledge of the complex molecular and cel- 
lular events during ABA signaling and of ABA-induced 
genes. Challenges and opportunities remain in exploring 
the complex interactions between ABA and defense in 
whole plants. 

Methods 

Plant material for RNA-Seq 

Seeds of tomato cv. Hongtaiyang 903 were planted and 
grown in plastic pots filled with organic loam in April, 
2012, and grown in a glasshouse in Chengdu, Sichuan 
Province, China. Sixty pots (four seedlings per pot) were 
used in this experiment. They were watered every other 
day. After 45 days when the plants had 5-7 leaves, the 
60 pots were divided into two groups. The plants in 30 
pots were sprayed with 400 mL of 7.58 i^mol L'^ ABA 



solution, and the plants in the other 30 pots were 
sprayed with the same volume of purified water as a 
control. Twenty-four hours later, the young third leaves 
of ten randomly-selected plants from both the ABA- 
treated and control groups were collected. The leaves of 
each group were combined, immediately snap-frozen, 
and stored in liquid nitrogen. These samples were la- 
beled aid (1 d after ABA-treatment) and eld (1 d, 
control). 

RNA extraction, cDNA library construction and lllumina 
deep sequencing 

Total RNA samples of aid and eld were prepared using 
Trizol reagent (Invitrogen, Carlsbad, CA, USA) and subse- 
quently used for mRNA purification and library construc- 
tion with the Truseq™ RNA Sample Prep Kit (lllumina, 
San Diego, CA, USA) following the manufacturer's in- 
structions. The two samples were sequenced on an lllu- 
mina HiSeq" 2000 (lllumina), generating 165,894,496 
reads. Each sample yielded more than 12 Gb of data. 
Sequencing was completed by the Shanghai Majorbio 
Bio-pharm Biotechnology Co. (Shanghai, China). 

Read trimming and optimization 

For each set of sequencing reads, the sequencing adapters 
were trimmed using SeqPrep (https://github.com/jstjohn/ 
SeqPrep), and then low quality bases (Solexa/IUumina 
quality score < 25) of the 3 ' ends were trimmed using in- 
house Perl scripts. The remaining high-quality reads were 
submitted for mapping analysis against the reference gen- 
ome sequences (ftp://ftp.solgenomics.net/tomato_genome/ 
annotation/ITAG2.3_release/) using Tophat [73]. The 
mapped reads were then assembled with Cufflink [29]. The 
assembled results and original genome annotations were 
merged and used for further annotation and differential- 
expression analysis. 

AH of these RNA-Seq reads were deposited in Sequence 
Read Archive database (http://www.ncbi.nlm.nih.gov/Traces/ 
sra/) under the Accession number of SRR926182 and 
SRR926185. 

Mapping reads to the reference genome and annotated 
genes 

Open reading frames (ORFs in all transcripts were pre- 
dicted using Trinity (http://trinityrnaseq.sourceforge.net/ 
analysis/extract_proteins_from_trinity_transcripts.html) 
[74] . Sequence-similarity Blast searches of these transcripts 
were conducted against the tomato genome reference, the 
NCBI NR protein database (http://www.ncbi.nlm.nih.gov/), 
the Gene Ontology (GO) database (http://www.geneontol- 
ogy.org/), the Search Tool for the Retrieval of Interacting 
Genes (STING) database (http://string-db.org/) [75], and 
the Kyoto Encyclopedia of Genes and Genomes (KEGG) 
database (http://www.genome.jp/kegg/) [76]. GO terms for 
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tomato transcripts were obtained using Blast2GO 
(v. 2.3.5) (http://www.blast2go.org/) with default parame- 
ters [30]. COG terms were obtained using Blastx 2.2.24+ 
in STRING 9.0. Metabolic pathway analysis was per- 
formed using Blastx/Blastp 2.2.24+ in KEGG (http://www. 
genome.jp/kegg/genes.html). 

Differential expression analysis 

The Tophat (http://tophat.cbcb.umd.edu/) and Cufflinks 
(http://cufflinks.cbcb.umd.edu/) programs provide FPKM 
(Fragments Per Kilobase of exon model per Million 
mapped fragments) values within a 95% confidence inter- 
val. Differential expression was analyzed and calculated ac- 
cording to the count values of each transcript in the two 
libraries using edgeR (the Empirical analysis of Digital 
Gene Expression in R) software [77]. "FDR < 0.05" and 
"|log2 fold-change (log2FC)| >1" were used as the thresh- 
olds for judging significant difference in transcript expres- 
sion. Transcripts with |log2FC| < 0.25 were assumed have 
no change in expression levels. 

qRT-PCR verification 

qRT-PCR was performed to verify the expression patterns 
revealed by the RNA-seq study. The purified RNA sam- 
ples were reverse-transcribed using the PrimeScript RT 
Reagent Kit with gDNA Eraser (Takara, Dalian, China) fol- 
lowing the manufacturer's protocol. Twenty-five tran- 
scripts were selected randomly for the qRT-PCR assay. 
Gene specific qRT-PCR primers (18-20 bp) (Additional 
file 1: Table S6) were designed using Premier 5.0 software 
(Premier Biosoft International, Palo Alto, CA). qRT-PCR 
was performed using SybrGreen qRT-PCR Master Mix 
(Ruian Biotechnologies, Shanghai, China) in an ABI 7500 
FAST Real-Time PCR System (Applied Biosystems, Foster 
City, CA, USA). PCR conditions were 2 min at 95°C, 
followed by 40 cycles of heating at 95°C for 10 s and an- 
nealing at 60°C for 40 s. Three replicates were performed, 
and the amplicons were used for melting curve analysis to 
check the amplification specificity. The relative expression 
level of each gene was calculated as 2'^^'^'^'^^ and the house- 
keeping gene glycemldehyde-3-phosphate dehydrogenase 
gene (Accession No. U93208) from S. lycopersicum was 
used to normalize the amount of template cDNA added in 
each reaction. 

Additional file 



Table S6. Comparison of expression patterns between RNA-seq expression 
and quantitative real-time PCR. 
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GAPDH: Glyceraldehyde-3-phosphate dehydrogenase; GLR: Glutaredoxin; 
GLU: Beta-l:3-glucanase; GO: Gene ontology; GPX: Glutathione peroxidase; 
GR: Glutathione reductase; GSH-AsA: glutathione-ascorbate; GST: Glutathione- 
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ORE: Open reading frame; PAL: Phenylalanine ammonia-lyase; 
PCR: Polymerase chain reaction; POD: Peroxidase; PP2Cs: Protein 
phosphatases which fall under the category of type 2C; PRO: Polyphenol 
oxidase; PRl; Pathogenesis-related proteins; PrxR: Peroxiredoxin; RCARs/PYRl/ 
PYL: Regulatory components of ABA receptor/pyrabactin resistance proteinl/ 
PYR-like protein; RNA: Ribonucleic acid; ROS: Reactive oxygen species; qRT- 
PCR: Quantitative real time - PCR; SA: Salicylic acid; sHSP: Small heat shock 
protein; SnRK2: The sucrose nonfermentingl -related protein kinase 2; 
SOD: Superoxide; STING: Search tool for the retrieval of interacting genes; 
TGA: TGACG/as-1: binding; Trx: Thioredoxin. 
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